Dynamics for partially coherent Bose-Einstein condensates in double wells 
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The dynamical properties of partially coherent Bose-Einstein condensates in double wells are 
investigated in three typical regimes. In the extreme Fock regime, the time evolution of the degree 
of coherence is shown to decay rapidly. In the Rabi regime, a relation between the amplitude of Rabi 
oscillation and the degree of coherence is obtained, which is expected to determine the degree of 
coherence by measuring the amplitude of Rabi oscillation. The study on the self-trapping phenomena 
in the Josephson regime exhibits that both the degree of coherence and the initial relative phase 
can affect the final particle distribution. 

PACS numbers: 03.75.Lm, 03.75.Hh 
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I. INTRODUCTION 

Quantum tunneling through a barrier, as a paradigm of 
quantum mechanics, has been observed in different sys- 
tems, such as two superconductors separated by a thin in- 
sulator [l| , two reservoirs of superfluid helium connected 
by nanoscopic apertures 0, and Bose-Einstein con- 
densates (BECs) trapped in double wells 3 ■ Comparing 
with the former two systems, BECs in double wells offer 
a versatile tool to study the quantum tunneling phenom- 
ena due to the fact that almost each parameter, such 
as the inter-well tunneling strength, the inter-particle in- 
teraction, and the energy bias between the two wells, 
can be tuned experimentally. The competition between 
those parameters makes BECs in double wells exhibit 
several fascinating phenomena, like Rabi oscillation, self- 
trapping, and Josephson oscillation, which have been ex- 
tensively studied 1, B S Hi- When a condensate is 
expected to be employed as a qubit, the first obstacle 
attempted to be avoided is the decoherence which was 
studied in experiment through the interference between 
BECs [9", ^Td\ . Theoretically, the effect of decoherence on 
the dynamics of BECs in double wells was recently dis- 
cussed with the help of single-particle density matrix [ll[ . 
The problem of decoherence in qubit measurement was 
recently investigated [l^l with a BEC in double wells. In 
current literature, the most of authors paid more atten- 
tion on the time evolution of the atom distribution but 
less attention on the relative phase between the conden- 
sates in two wells. Meanwhile, the systems were mostly 
assumed to be completely coherent at the initial time. 
Since both relative phase and degree of coherence are 
believed to affect the dynamical properties of BECs in 
double wells, it is worthwhile to study the coherence dy- 
namics with attention to the relative phase and the de- 
gree of coherence. 

In this paper, we study a BEC system in double wells 
with different degrees of coherence. In mean-field ap- 
proximation, we study the dynamical properties of the 
system in different regimes with the help of the single- 
particle density matrix and show that, in comparison to 
the completely coherent case, partially coherent BECs 



in double wells can exhibit richer physics. In the next 
section, we model the partially coherent system and give 
the mean-field dynamical equation for the elements of 
single-particle density matrix. In Sec. IIIIl we study the 
evolution of the degree of coherence and discuss the dy- 
namical property of the partially coherent system in Fock 
regime. In Sec. IIVI we study the system in Rabi regime 
and discuss the influence of degree of coherence on the 
Rabi oscillation. In Sec. |Vl we study self-trapping phe- 
nomenon for the partially coherent system. Then brief 
summary and discussion are given in Sec. IVII 



II. MODELING PARTIALLY COHERENT 
SYSTEMS 

We consider a Bose-Einstein condensate confined in a 
double-well potential, where atoms can tunnel between 
the two wells. The Hamiltonian of such a system in the 
second-quantization form is given by, 



H = -In, 



(1) 



where the bosonic operators ajj and — 1,2) creates 

and annihilates an atom in the /ith well, respectively; and 
fi^ = ajj^a^ is the particle number operator. Here the pa- 
rameter 7 denotes the energy bias between the two wells, 
T is the inter- well tunneling strength, and U is the in- 
teraction strength between atoms. The Hamiltonian ([T]) 
can describe not only the BEC in double wells but also 
that in two hyperfine states [H, This model was 

conventionally studied in a mean-field approach by re- 
placing the expectation values of annihilators in two dif- 
ferent wells with two complex numbers oi and 02, respec- 
tively 15, 16] . In those works, the system is essentially 
assumed to be of complete coherence, i.e., the conden- 
sate stays in a completely coherent superposition state 
I V'coh) — -/^(aia^ + a2al)^ \ vac). Whereas, a realistic 
system of condensates may be not always in complete co- 
herence for various situations. It is therefore worthwhile 
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to study the dynamical properties of partially coherent 
BECs in double wells. 

As we know, it is convenient to introduce the 
single-particle density matrix p with entities p^^{t) = 
( djj^{t)a^{t) ) where the expectation value is taken for the 
initial state of the system. Clearly pu and P22 represent 
the population in the first and in the second well, respec- 
tively. With the help of Heisenberg equation of motion 
for operators, one can derive the dynamical equations for 
the elements of the aforementioned density matrix by the 
mean- field approach in the semiclassical limit. 



■ dPii 

dt 

■ dpi2 
'' dt 

■ dp2i 
' dt 



; dp22 
* dt 



'T{Pi2 - P21). 



-7Pl2 + T{P22 - Pll) - UN{p^-^ - P22)Pl2: 
7P21 - T{P22 - Pll) + UN{p^-^ - P22)P21. 



(2) 



where h is set to unit and p^^ = p^^/iV. The conservation 
of particle number requires that Pn ^- P22 — 1- 

As the 2x2 density matrix can be expanded as p — 
{I + P ■ a)/2 with (t's the Pauli matrices and P a vector 
inside the so-called Bloch sphere. Since |P| = 1 refers to 
a pure state (completely coherent superposition) while 
|P| < 1 refers to a mixed state (partial coherence), one 
can measure the degree of coherence by jPp = 77. Thus 
a natural definition of the degree of coherence is given 
by 



7] = 2Trp2 - 1, 



(3) 



which is an important quantity that affects the dynam- 
ical features, such as Rabi oscillation, self-trapping, etc. 
With the help of the mean-field dynamical Eqs. ^ for 
the single-particle density matrix, the dynamics of the 
system can be investigated. We know that the ratio of 
the interaction strength U to the tunneling strength T de- 
termines three distinct regimes, namely, the Rabi regime, 
UN/T < 1, the Josephson regime, 1 < UN/T < N^, 
and the Fock regime, iV^ ^ UN/T. The system man- 
ifests different dynamical features in different regimes, 
which will be given in the following sections. 



III. EVOLUTION OF THE DEGREE OF 
COHERENCE IN THE FOCK REGIME 

To study the evolution of the degree of coherence, we 
can conveniently introduce the pseudospin operators. 



1 



Jz = ^{a\di - dla^), 
Jx = ]^{d\d2 + a2«i)' 



J u 



-{d\a2 - d\a{), 



which obey commutation relations for the angular mo- 
menta [Jj, JjJ\ = ie.^^Ji, and fulfil 



N /N 



2 V 2 



1 



for systems of N bosons. In terms of these pseudospin 
operators, the Hamiltonian ([T]) can be written as, 



2T.h 



(4) 



This implies that the dynamical properties are deter- 
mined only by the direction of the pseudospin J since 
its magnitude is fixed on N/2. 

We know that the eigenvalue of Jz is (ni — n2)/2, which 
refers to the population imbalance between the two wells. 
For convenience, let us introduce the canonically conju- 
gated operator ip of Jz to characterize the relative phase, 
which can be regarded as the angle of J in the x-y plane in 
the angular-momentum picture. According to Ref. p7| . 
the operator ip can be defined through E = exp i(p, where 

E = [(iV/2 - Jz){N/2 + Jz + + iJy). (5) 

Such a definition satisfies [Jz , E] = exp iip which is con- 
sistent with the condition [Jz,<^] = —i, so that Jz and (p 
are canonically conjugated to each other. Then, in terms 
of Jz and p, equation (|4]) can be rewritten as 



H = UJt+ jJz - TN\I 1 - ^ cos p. (6) 



Uz 



In the Fock regime, the tunneling strength is much 
smaller than the interaction one between atoms, i.e., 
U/{NT) > 1, hence the last term in Eq. ([6]) can be 
neglected. Such a condition can be satisfied in exper- 
iment through increasing the distance between the two 
wells or enhancing the height of the potential barrier sep- 
arating the two wells. Then the dynamical equations in 
this regime become. 



dJz 
dt 



0, 



dp 
It 



2UJz + 7, 



which means that the difference in particle numbers 
between the two wells is fixed but the relative phase 
p = 2UJzt + "ft evolves with time. According to Eq. ([3]) 
and the definition of pseudospin operators, the degree of 
coherence can be written as 



4 - 



d\d2)?) 



From the definition of p, one can also get 



d\a2 



Since the expectation value of Jz is fixed in the extreme 
Fock limit T = 0, the evolution of the degree of coherence 
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only depends on that of (p. For example, if the particle 
numbers in the two wells are the same at the initial time, 
the expectation value of Jz will be always fixed on zero in 
the extreme Fock limit. In this case, the degree of coher- 
ence becomes = 4|( d|a2 )P/-/V^ « |( expi<^ )p, which 
reflects that the evolution of the degree of coherence is 
determined by the interaction strength f7, the detuning 
7 and the initial state of the system. 

As we know, the Fock bases 1 1) composing the (A^ + 1)- 
dimensional Hilbert space of the system are 



N N 



J _ N N 



1, 



N 



(7) 



where / corresponds to the quantum number of Jz de- 
noting the half of the difference in the particle numbers 
between the two wells. Since any initial state can be 
written as 

l=N/2 

i*(o))= Yl ^'10, 

the degree of coherence can be obtained as long as tpi is 
given. Taking ipi = e~^''~^^ / (^1^/2)^^^ as an example 
and replacing the sum over / by an integral, one can 
obtain the degree of coherence. If keeping the lowest 
order term in the Taylor expansion in the calculation of 
the expectation value of (1 — -J^- J|)^/^ expi(p, we have 



7V2 



(8) 



This expression is valid when the width of the Gaus- 
sian distribution is much smaller than TV, i.e., Iq N. 
Clearly, the degree of coherence decays with time as long 
as there is interaction between atoms. Such result tells 
us that one can prepare systems with different degrees of 
coherence through changing the evolution time t in the 
extreme Fock regime. From Eq. (|8]), we can also find that 
the larger the is, the faster the degree of coherence will 
decay. Then one can suppress the decay through decreas- 
ing the value of Iq. Note that in the above calculation of 
Tjit), we used = 2UJzt + jt which is valid only for the 
extreme Fock regime. In the other regimes, one needs to 
solve Eq. ^ either analytically or numerically without 
any assumption. 



IV. THE INFLUENCE OF DEGREE OF 
COHERENCE ON THE RABI OSCILLATION 

Rabi oscillation is an important phenomenon reflecting 
the coherent property of a system, which has been dis- 
cussed recently 8_, 18, 19]. In Ref. 8] the initial state was 
assumed to be a pure state, i.e., pu — 1, implying r] — 1. 
Whereas, according to our formulation in Sec. Illli one 
can prepare states with different degrees of coherence, 
which makes Rabi oscillation worthwhile to be studied 
from a new angle of view. Now we study the dynamics of 



Rabi oscillation in terms of the density matrix and give 
the relation between the amplitude of the Rabi oscillation 
and the degree of coherence for BECs in double wells. 




FIG. 1: (Color online) The Rabi oscillation of the particle 
distribution with different degrees of coherence for (a) U /T = 
and for (b) UNjT = 4. The parameters are 7/T = 0.5 and 
r\ = 0.5 (dot line), 1 (solid line). 

In the extreme Rabi limit, IJ = 0, equations 
are reduced to Bloch equations that can be solved an- 
alytically. For example, we consider a initial state de- 
scribed by the density matrix Pii(O) — ^22(0) = 1/2 
and Pi2(0) = P2i(0) — c, we can obtain the solution of 
Eqs. @, 



P22 
P12 
P21 



1 2c-/T 2c-fT 

2 ^ ~W ^ ~W 
1 2c^T 2cjT 

2^~f^ 



cos (r^i), 

2 COS (fit). 



9 ? 

c - + COS iSlt) -f i-^ sin (Sit), 



C7 

02" 



C7 
1^ 



COS (rjt) - i^sin(m), (9) 



where ^ = {^T'^+^i^ f^'^ , and c is a real number related to 
the degree of coherence, i.e., c = ^/ri/2. From Fig. [l] we 
can see that the particle numbers in the two wells both 
periodically oscillate with time. The oscillation ampli- 
tude c-^T 1^ diminishes with the decrease in the degree 
of coherence. Here the relative phase between the con- 
densates in two wells is taken as zero, i.e., pi2(0) is real 
at the initial time. If the degree of coherence is fixed but 
the relative phase is not zero, saying pi2(0) = cexp(i0) 
and P2i(0) = cexp(— i0), we can also solve Eqs. ^ an- 
alytically (the expression of Pij(t) is omitted for saving 
space). One can find that the amplitude of Rabi oscil- 
lation is affected by the initial values of the degree of 
coherence as well as the relative phase. For the sake of 
comparison, we plot a numerical result of the Rabi oscil- 
lation in the presence of interaction {U ^ 0) in Fig.[T](b). 
Comparing the two panels in Fig. [l] we can see that the 
interaction between atoms will make both the amplitude 
and the period of the Rabi oscillation decrease for the 
same initial states. 

The above discussions tell us that one can obtain the 
degree of coherence through measuring the amplitude of 
the Rabi oscillation in experiment. For example, like in 
the experiment f 20j , prepare a BEC in one hyperfine state 
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and transfer one half of atoms into the other hyperfine 
state through a two-photon pulse; turn off the pulse to 
allow the system to evolve freely without the inter-state 
tunneling (i.e., T = 0) until t — tq. Then turn on a pulse 
which makes the atoms tunnel between the two states and 
measure the amplitude of the Rabi oscillation. Because 
the system is in the extreme Fock regime when t < tq, 
based on the calculation given in Sec. IIIIl we can have 

Pii = P22 = 1/2 and p^^^ 2 '^^P^*'^'^") ^^^^"'o^^'^o 
at the time tq. Such a state can be actually prepared as 
the initial state of the subsequent Rabi oscillation proce- 
dure. Note that the expression of pi2 can be also rewrit- 
ten as pi2 = Y^exp(i7To)/2 considering the definition of 
r]. Since the relative phase jtq is determined, the am- 
plitude of Rabi oscillation only depends on the degree of 
coherence. The analytical result of Eqs. ^ can not be 
obtained due to the existence of interaction terms in the 
above example, so one can solve it numerically. The de- 
gree of coherence can be determined through fitting the 
experimental Rabi oscillation profiles with the theoretical 
ones as the degree of coherence is evolved. 



V. SELF-TRAPPING FOR THE PARTIALLY 
COHERENT SYSTEM 

As we know, the most prominent feature of atomic tun- 
neling between two wells is the nonlinear dynamics aris- 
ing from the interaction of atoms. Whereas, equations ([2|) 
can not be analytically solved once the interaction terms 
are taken into account. In this case, we solve Eqs. ([2]) 
numerically. In the numerical calculation, we adopt a 
linearly time-dependent energy bias "f = at where a is 
a constant characterizing the rate of the change of the 
energy bias 7. The initial values are pn = P22 = 1/2 

and Pi2 — —y/rjexpi(j3, where (j) is the phase difference 

between the condensates in the two wells at initial time. 
Our results manifest that the initial relative phase and 
the degree of coherence affect the dynamical properties 
of self-trapping. 

In Fig. [21 we plot the time-evolution of the population 
Pii for different initial values of the degree of coherence 
and relative phase. From this figure, we can see that the 
system exhibits the phenomenon of self-trapping if the 
system is coherent, i.e., 77 ^ 0, but the atoms do not tun- 
nel between the two wells for the incoherent system, i.e., 
T] = 0. We can also find that most of the atoms favorite 
to stay in the right well (i.e., P22 > Pu) at the final time 
for (j> = 7r/3, which is contrast to the case for = tt. The 
difference between the two panels of Fig. [2] implies that 
the initial relative phase can affect the phenomenon of 
self-trapping. It depends on the energy bias and the ini- 
tial relative phase that in which well the atoms prefer to 
stay at the final time. When the initial relative phase is 
fixed, figure [2] shows that the larger the degree of coher- 
ence is, the larger the population difference between the 
two wells at the final time will be. This is confirmed by 




' t(1/T) » 



FIG. 2: (Color online) Time-evolution of the population prob- 
ability for different degrees of coherence and initial relative 
phases ((^ = 7r/3 for the top panel and <^ = tt for the bottom 
panel). The other parameters are UN/T — 4 and a/T^ — 5. 
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FIG. 3: (Color online) The dependence of the final population 
probability on the degree of coherence. The parameters are 
the same as in Fig. (2] 



Fig. [3] where the final population probability versus the 
degree of coherence is plotted for three different values 
of (j). In Fig. m we plot the relation between the final 
population probability and the initial relative phase for 
different degrees of coherence. This figure confirms that 
the initial relative phase can affect the final distribution 
of atoms. 

The emergence of the self-tapping phenomena is de- 
pendent on the initial state and the systems param- 
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bias, the system can not exhibit the self-trapping phe- 
nomena due to the critical interaction strength Uc oo 
for z(0) = 0, 0(0) ^ ttIH. For the initial states with 
population imbalance {i.e., z(0) ^ 0), we plot the nu- 
merical solutions of Eqs. ([2]) for a symmetric double- well 
potential with different degrees of coherence in Fig. [5l 
Such a figure shows that the degree of coherence affects 
the critical interaction strength Uc above which the sys- 
tem can exhibit the self-trapping phenomena in a sym- 
metric double-well potential. The system investigated in 
Ref. Q corresponds to the case when the degree of coher- 
ence = 1 in our discussion (see the dot line in Fig. [5]). 



FIG. 4: (Color online) The dependence of the final population 
probability on the initial relative phase (p. The parameters are 
the same as in Fig. (2] 



VI. SUMMARY AND DISCUSSION 
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FIG. 5; (Color online) Time-evolution of population imbal- 
ance for a symmetric double-well potential with different de- 
grees of coherence. The initial states are pn — 0.8, p22 = 0.2, 
and pi2 = 0.4 (dot line), pi2 = 0.2 (solid line). The parame- 
ters are (a) UN/2T = 1, (b) UN/2T = 8, (c) ;7iV/2r = 9.99, 
and (d) UN/2T = 11. 



eters, e.g., the interaction strength U, the inter-well 
tunneling strength, and the energy bias 7 between two 
wells. Such a phenomenon was investigated in a sym- 
metric double- well potential i.e., 7 = 0, and also 
in a double-well potential with a periodic modulation, 
i.e., 7 cx sinwi [2l|. According to Ref. [5|, we can 
find that the self-trapping phenomena occur only when 
the interaction strength is larger than the critical value 
Uc oc [^1 - z(O)2cos0(O) + l]T/[Nz{0)^] in the case of 
7 = 0, where z(0) refers to the initial population differ- 
ence pii(O) — P22(0) and 0(0) the initial relative phase 
between the condensates in the two wells. Therefore, for 
the initial states we considered afore, if there is no energy 



In the above, we considered Bose-Einstein condensates 
in double wells with different degrees of coherence. With 
the help of single-particle density matrix, we studied the 
dynamical properties of partially coherent Bosc-Einstcin 
condensates in double wells and showed that the degree 
of coherence is a useful parameter that affects the dy- 
namical features. We investigated the system in different 
regimes and found that the degree of coherence can affect 
the dynamical properties of the system significantly. In 
the Fock regime, we mainly studied the time evolution 
of the degree of coherence by introducing the pseudospin 
operators and showed that the degree of coherence de- 
cays in exponential form of the square of time. In the 
Rabi regime, we studied the effects of the initial rela- 
tive phase and degree of coherence 77 on Rabi oscilla- 
tion and showed that the amplitude of Rabi oscillation 
is in proportional to the square root of rj in the case 
of (f> — 0. According to the relevant result, the degree 
of coherence is expected to be determined through mea- 
suring the amplitude of Rabi oscillation. Because the 
existence of nonlinear terms of interaction makes the dy- 
namical equations in the Josephson regime not solvable 
analytically, we solved those equations numerically and 
found that the self-trapping phenomenon also exists for 
partially coherent BEC systems in double wells. A more 
fascinating feature is that which well the particles will 
stay in at the final time largely depends on the initial 
relative phase. This result is different from that in the 
previous works where the particles are in the same 
well at the initial time such that the initial relative phase 
between BECs in two wells does not affect the dynamical 
property of systems explicitly. 
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